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1. 


Introduction 


The  potential  instability  of  disposal  orbits  of  the  Global  Positioning  System  (GPS)  was  discovered  by  Chao 
(2000).  The  instability  manifests  itself  as  significant  growth  in  orbital  eccentricity  over  a  timeframe  of 
decades.  The  cause  of  this  long-term  eccentricity  growth  is  a  dynamical  resonance  condition  resulting  from 
the  combined  gravitational  pull  of  the  Sun,  Moon,  and  the  non-spherical  gravity  field  of  the  Earth.  Long-term 
eccentricity  growth  of  disposal  orbits  will  lead  to  penetration  of  the  constellation  by  previously  disposed 
vehicles,  thereby  posing  a  collision  risk  to  the  operational  constellation  members. 

Chao  and  Gick(2002)  also  showed  that  disposal  orbits  for  Glonass  and  Galileo  can  be  similarly  unstable. 
Hence,  long-term  disposal  orbit  eccentricity  growth  affects  all  of  the  existing  and  planned  satellite 
constellations  in  medium  Earth  orbit  (MEO). 

Gick  and  Chao  (2001)  showed  that  the  amount  of  eccentricity  growth  depends  on  the  initial  elements  of  the 
disposal  orbit.  In  particular,  minimization  of  the  disposal  orbit  initial  eccentricity,  e0,  and  proper  selection  of 
the  initial  argument  of  perigee,  cd0,  can  suppress  eccentricity  growth,  and  hence  constellation  penetration, 
over  a  time  period  of  up  to  200  years.  The  initial  right  ascension  of  ascending  node  (RAAN  =  Q0)  and 
inclination,  i0,  also  strongly  influence  long-term  eccentricity  growth  (Chao  and  Gick,  2002),  but  these  orbital 
parameters  are  not  easily  modified  during  the  disposal  process  due  to  the  excessive  amount  of  required  AV. 

Jenkin  and  Gick  (2001,  2002)  determined  the  long-term  collision  risk  posed  to  the  operational  GPS 
constellation  by  disposed  vehicles  placed  in  orbits  with  perigees  500  km  and  832  km  above  semi- 
synchronous  orbit.  In  these  cases,  the  disposal  strategy  that  was  modeled  involved  achieving  a  low  initial 
eccentricity  of  0.005.  The  initial  argument  of  perigee  was  not  selected  to  minimize  eccentricity  growth,  but 
instead  to  yield  average  eccentricity  growth  after  200  years.  The  intent  was  to  model  the  average  collision 
risk  that  would  result  if  co0  were  allowed  to  take  on  random  values. 

Jenkin  and  Gick  (2003)  also  determined  the  long-term  collision  risk  posed  to  the  operational  GPS 
constellation  by  disposed  upper  stages  that  cannot  raise  their  disposal  orbit  perigees  outside  of  the 
constellation.  The  disposal  orbits  of  these  upper  stages  have  relatively  high  apogees  above  the  constellation, 
and  therefore  significant  initial  eccentricity.  As  a  result,  the  upper  stage  disposal  orbits  undergo  more  rapid 
eccentricity  growth  than  would  occur  if  the  vehicles  were  disposed  with  perigee  500  km  above  the 
constellation,  with  low  e0  and  random  co0.  The  study  showed  that  the  early  collision  risk  posed  to  the 
constellation  by  the  upper  stage  disposal  orbits  with  perigee  left  in  the  constellation  will  be  higher  than  if  the 
vehicles  were  disposed  with  perigee  500  km  above  the  constellation  (with  low  e0  and  random  co0)-  However, 
eventually  the  collision  risk  posed  by  the  case  of  vehicle  disposal  at  500  km  would  catch  up  and  exceed  the 
collision  risk  posed  by  the  upper  stage  disposal  orbits.  It  was  evident  that  the  faster  eccentricity  growth  of  the 
upper  stage  disposal  orbits  results  in  a  dilution  of  the  long-term  collision  risk. 

The  purpose  of  the  study  presented  here  is  to  investigate  the  possibility  of  diluting  MEO  disposal  orbit 
collision  risk  by  selection  of  e0  and  co0  in  order  to  control  long-term  eccentricity  growth.  The  Galileo 
constellation  was  selected  as  an  example. 
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2. 


Overall  Methodology 


The  study  analysis  flow  was  as  follows.  The  long-term  evolution  of  the  expected  range  of  constellation  orbits 
was  simulated  using  a  mean  element  propagator  of  high  accuracy  to  account  for  eccentricity  growth.  From 
the  resulting  data,  a  statistical  spatial  distribution  model  of  the  operational  constellation  was  generated.  Next, 
for  each  disposal  strategy  considered,  the  entire  expected  range  of  disposal  orbit  initial  conditions  was  swept 
out,  and  the  resulting  disposal  orbits  were  propagated  over  200  years.  The  propagated  disposal  orbit  histories 
and  the  constellation  spatial  density  model  were  then  used  to  generate  an  ensemble  of  time  profiles  of 
cumulative  long-term  collision  probability,  one  for  each  disposal  orbit  initial  condition.  The  total  collision 
risk  accounting  for  on-going  disposal  of  satellites  over  time  was  determined  by  randomly  selecting  collision 
probability  time  profiles  from  the  ensemble  that  was  generated,  time-shifting  them  to  account  for  future 
disposal  epochs,  and  then  summing  them  together.  The  entire  disposal  sequence  over  200  years  is  repeated 
1000  times  in  Monte  Carlo  fashion.  The  overall  methodology  is  depicted  in  Figure  1. 


Figure  1.  Overall  study  methodology. 
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3. 


Long-Term  Orbit  Propagation 


The  Aerospace  Corporation  tool  MEANPROP  was  used  to  perform  the  long-term  propagation  of  the 
constellation  and  disposal  orbits.  MEANPROP  is  a  mean  orbit  element  control  simulation  that  uses  the  Semi- 
Analytic  Orbit  Propagator  (SAOP)  to  perform  long-term  propagation.  SAOP  is  a  program  developed  by  the 
Charles  Stark  Draper  Laboratory  that  has  undergone  extensive  validation  (McLain,  1977).  In  this  study,  the 
force  model  included  Sun-Moon  gravity,  a  12  x  12  WGS84  Earth  gravity  field,  solar  radiation  pressure,  and 
atmospheric  drag.  For  modeling  the  effect  of  solar  radiation  pressure,  vehicle  design  data  on  the  GSTB-V2-B 
experimental  Galileo  satellite  was  taken  from  Space  Daily  (2004).  As  a  result,  the  vehicle  mass  was  assumed 
to  be  523  kg,  projected  average  cross-sectional  area  was  assumed  to  be  4.53  m2  (obtained  by  prorating  GPS 
area  by  the  Galileo  power  level),  and  the  reflectivity  coefficient  cr  was  assumed  to  be  1.3.  For  modeling  the 
effect  of  atmospheric  drag,  the  MSIS-90  atmosphere  model  was  used,  the  solar  flux  parameter  FI 0.7  was  set 
to  a  constant  value  of  140,  the  geomagnetic  index  Ap  was  set  to  a  constant  value  of  10,  and  the  drag 
coefficient  was  assumed  to  be  2.0. 
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4. 


Collision  Risk  Analysis 


A  density-based  method  was  used  to  compute  collision  risk  in  this  study.  This  method  is  formulated  from  the 
perspective  of  the  primary  satellite  as  it  flies  through  a  field  of  secondary  objects.  The  following  formulation 
was  used: 


Nc  (t)  =  Acc  J  p(rs  (r))  v  dr  (i} 

*0 


Pc(t)  =  l-e~Nc(t)  (2) 


where  t  is  the  current  time  point,  t0  is  the  propagation  start  epoch,  pc(t)  is  the  cumulative  collision  probability 
at  time  t,  Nc(t)  is  the  average  number  of  collisions  at  time  t,  rs  is  the  position  of  the  primary  satellite,  p(rs)  is 
spatial  density  (averaged  over  latitude  and  time)  of  the  secondary  objects  at  the  position  of  the  primary 
satellite,  v  is  the  average  relative  velocity  between  the  primary  and  secondary  satellites,  and  Acc  is  the 
average  collision  cross-sectional  area  between  the  primary  and  secondary  satellites. 

The  average  collision  cross-sectional  area  Acc  was  set  to  unity,  thereby  normalizing  the  collision  probability. 
This  was  done  because  little  information  was  available  in  the  open  literature  on  the  Galileo  spacecraft  design. 
Since  the  focus  of  this  study  is  a  comparison  of  the  effectiveness  of  different  disposal  strategies  on  collision 
probability,  the  absolute  level  of  the  collision  probability  is  not  necessary.  It  is  sufficient  to  compare  collision 
probabilities  normalized  by  cross-sectional  area.  It  is  a  simple  matter  to  obtain  the  absolute  collision 
probability  level  by  scaling  the  normalized  probabilities  with  actual  collision  cross-sections  determined  from 
design  information. 

In  order  to  accelerate  computer  processing,  the  integration  in  Eq.  1  was  implemented  using  the  method  of 
averaging.  This  involves  separating  the  problem  into  a  short  time-scale  problem  and  a  long  time-scale 
problem.  The  short  time-scale  problem  consists  of  determining  the  average  number  of  collisions  over  a  single 
orbit  revolution  of  the  primary  satellite  for  an  instance  of  mean  orbital  elements  that  have  been  propagated  to 
a  given  epoch  by  MEANPROP.  The  average  number  of  collisions  per  revolution  is  then  divided  by  the  orbit 
period  to  yield  the  average  rate  of  growth  of  average  number  of  collisions  at  the  given  mean  orbital  element 
epoch. 


) =  y  A,:c  { p^s  (Zs  ^vd  Ts  (3) 

where  Tp  is  the  orbital  period  of  the  primary  satellite.  The  long  time-scale  problem  then  consists  of 
integrating  Eq.  3. 


A'C(<I)  =  J 

*0 


{rL)drL 


(4) 
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A  computer  program  was  used  that  evaluates  Eq.  3  at  each  time  point  in  the  mean  orbit  element  histories 
(i.e.,  at  each  long  time-scale  point).  It  then  evaluates  the  integral  in  Eq.  4  (and  subsequently  Eq.  2)  via  a 
discrete  trapezoidal  algorithm  that  processes  the  results  of  Eq.  3  at  those  time  points.  This  program  was  used 
to  generate  profiles  of  cumulative  collision  probability  vs.  time  for  each  instance  of  initial  disposal  orbit 
conditions. 

The  density-based  method  used  in  this  study  was  selected  because  it  is  computationally  fast.  However,  it  is 
an  approximate  technique  that  is  less  accurate  than  the  miss  distance-based  method  that  has  been  used  by  the 
authors  (Jenkin  and  Gick,  2001,  2002).  It  does  not  account  for  the  correlations  that  exist  between  the  disposal 
orbit  and  the  constellation  satellite  orbits.  In  particular,  there  are  correlations  between  the  two  orbit  classes  in 
Q  and  oo  because  they  evolve  slowly,  especially  relative  to  each  other.  However,  the  values  of  Q  for  the 
constellation  planes  are  uniformly  spaced  as  a  whole  (albeit  deterministically  and  not  randomly),  and  the 
initial  values  of  co  of  the  constellation  slots  are  randomly  distributed  throughout  the  constellation.  As  a  result, 
the  accuracy  loss  of  the  density-based  method  is  not  severe. 

To  demonstrate  this  fact,  a  numerical  comparison  between  the  two  methods  was  performed  for  a  GPS 
disposal  case  taken  from  Jenkin  and  Gick  (2002).  In  this  case,  six  vehicles  are  disposed  at  500  km,  one  in 
each  constellation  plane.  These  vehicles  subsequently  penetrate  the  GPS  constellation  due  to  orbital 
eccentricity  growth.  Using  each  method,  the  collision  probability  time  profiles  over  200  years  were 
computed  and  then  averaged  together  into  a  single  profile.  The  averaged  profiles  resulting  from  the  two 
methods  are  presented  in  Fig.  2.  This  plot  shows  that,  for  the  case  considered,  the  density-based  method 
generated  a  collision  probability  profile  that  is  higher  than  that  computed  by  the  miss  distance-based  method 
by  a  factor  of  at  most  1.4.  Given  other  uncertainties  in  the  problem  (constellation  configuration,  satellite 
replacement  rates,  disposal  orbit  initial  conditions,  collision  radii),  this  deviation  is  considered  to  be 
acceptable  for  this  study. 


3.0E-06 


2.5E-06 
S  2.0E-06 

-Q 

O 

cl  1.5E-06 

c 

o 

§  1.0E-06 

o 

O 

5.0E-07 


0.0E+00 

0  20  40  60  80  100  120  140  160  180  200 


Time  (years) 

Figure  2.  Comparison  of  methods  for  computing  collision  risk.  Results  are  for  a  GPS  disposal  case  taken  from  Jenkin 
and  Gick  (2002). 
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5. 


Constellation  Model 


In  order  to  account  for  the  effect  of  eccentricity  growth  on  the  orbital  configuration  of  the  constellation,  the 
mean  elements  of  constellation  orbits  were  propagated  over  an  assumed  lifespan  of  the  satellites.  It  was 
assumed  that  satellite  station  keeping  maneuvers  are  carried  out  in  a  way  so  that  they  do  not  directly  change 
eccentricity.  This  assumption  is  taken  from  the  practice  for  GPS  (Jenkin  and  Gick,  2001).  The  effect  of 
station  keeping  on  in-track  position  was  not  modeled.  Therefore,  the  longitude  of  ascending  node  was 
allowed  to  drift  in  the  propagation  runs.  It  is  not  necessary  to  model  the  effect  of  station  keeping  on  the  in¬ 
track  position  because  the  in-track  position  has  no  effect  on  long-term  collision  probability.  The  formulation 
of  collision  probability  used  in  this  study  (described  below)  is  based  on  the  assumption  that  the  in-track 
position  (i.e.,  argument  of  latitude)  of  the  disposed  vehicle  relative  to  the  satellite  is  continually  randomized 
over  360  deg.  This  assumption  is  valid  due  to  the  difference  between  the  semimajor  axes  of  the  disposed 
vehicle  and  the  constellation  satellites.  This  difference  in  semimajor  axes  causes  a  large  enough  difference  in 
orbital  mean  motions  to  produce  continual  re-circulation  of  the  disposed  vehicle  relative  to  the  satellite  over 
360  deg  with  a  repeat  period  of  days  to  weeks  (depending  on  the  disposal  orbit  semimajor  axis). 

The  propagation  start  epoch  is  assumed  to  be  1  January  2009.  The  initial  semimajor  axis  of  constellation 
orbits  was  assumed  to  be  29994  km,  and  inclination  was  assumed  to  be  56  deg  (ESA,  2002).  The  orbital 
eccentricity  was  assumed  to  have  an  initial  value  of  0.008,  which  is  taken  from  GPS  practice  (Jenkin  and 
Gick,  2002). 

The  constellation  orbits  will  undergo  nodal  regression  due  to  the  J2  oblateness  term  of  the  Earth’s 
gravitational  field.  The  nodal  regression  rate  for  the  Galileo  constellation  is  9.025  deg/year.  Together,  all 
three  constellation  planes  will  sweep  out  the  entire  360  deg  range  of  right  ascension  of  ascending  node,  Q,  in 
13.3  years.  In  addition,  the  argument  of  perigee,  oo,  may  take  on  any  value  over  its  360  deg  range.  Therefore, 
to  generate  the  constellation  model,  the  initial  values  Q0  and  co0  were  varied  over  the  range  [0,360)  deg  in  10 
deg  increments.  (The  notation  [0,360)  deg  denotes  the  entire  range  from  0  deg  to  360  deg,  including  0  deg 
but  excluding  360  deg  because  it  is  redundant  with  0  deg.)  Orbital  mean  element  vs.  time  profiles  for  each 
(Q0,  co0)  pair  were  generated  by  propagating  over  a  time  period  of  15  years  using  MEANPROP.  This  time 
period  is  based  on  the  assumption  that  the  maximum  vehicle  life  is  15  years. 

Figure  3  shows  the  resulting  evolution  of  apogee  and  perigee  altitude  (relative  to  the  constellation  reference 
orbit)  for  a  sample  set  of  constellation  orbits  with  Q0  =  180  deg.  Each  curve  corresponds  to  a  specific  value 
of  co0. 
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Figure  3.  Apogee  and  perigee  altitude  profiles  of  Galileo  constellation  orbits  (shown  relative  to  the  constellation 
reference  orbit).  Curves  are  shown  for  Q0  =  180  deg  and  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg 
increments. 


A  statistical  model  of  the  spatial  distribution  of  the  constellation  satellites  was  generated  to  support  the 
collision  risk  analysis  in  this  study.  More  specifically,  the  spatial  number  density  as  a  function  of  altitude  was 
computed.  This  distribution  was  determined  as  follows.  The  spatial  density  distribution  over  altitude  at  each 
time  point  in  the  MEANPROP  long-term  propagation  files  was  computed  using  the  formulation  of  satellite 
altitude  distribution  derived  by  Dennis  (1972).  This  formulation  is  equivalent  to  the  one  derived  by  Kessler 
(1981).  The  resulting  spatial  density  distributions  for  each  time  point  are  intrinsically  averaged  over  latitude. 
These  spatial  density  distributions  were  then  averaged  across  all  the  propagation  time  points  and  (Q0,  ooo) 
pairs,  and  then  prorated  by  the  number  of  constellation  members,  which  is  assumed  to  be  30  (ESA,  2002). 
The  resulting  constellation  spatial  density  distribution  assumes  that  the  Q  and  co  for  the  orbits  are  randomly 
and  uniformly  distributed  over  360  deg.  This  model  therefore  is  a  statistically  stationary  (i.e.,  static) 
representation  of  the  constellation  that  accounts  for  replacement  of  satellites  in  a  random  order  at  an  average 
rate  of  2  satellites  per  year  (  =  30  constellation  satellites  15  year  satellite  life). 

Figure  4  shows  the  resulting  average  spatial  density  distribution.  It  can  be  seen  from  this  figure  that  altitude 
boundaries  at  +/-  500  km  (relative  to  the  constellation  reference  orbit)  effectively  clear  the  Galileo 
constellation  most  of  the  time.  A  disposal  keepout  zone  of  +/-  500  km  is  a  useful  reference  because  it  is  the 
keepout  zone  that  is  specified  for  the  semi-synchronous  orbit  by  the  U.S.  Government  Debris  Mitigation 
Standard  Practice. 
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Altitude  relative  to  Galileo  ref.  orbit  (km) 

Figure  4.  Variation  of  spatial  number  density  with  altitude  (relative  to  the  constellation  reference  orbit)  of  the  Galileo 
constellation  model  used  in  this  study.  The  spatial  density  is  averaged  over  latitude. 
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6.  Concept  Of  Collision  Risk  Dilution 


In  order  to  understand  the  concept  of  collision  risk  dilution,  the  dependence  on  orbital  eccentricity  of 
probability  of  collision  posed  by  a  sample  disposed  Galileo  satellite  was  determined.  The  disposed  satellite  is 
assigned  a  semimajor  axis  a  =  30647  km.  When  the  eccentricity  is  0.005,  the  perigee  is  500  km  above  the 
Galileo  reference  circular  orbit. 

The  collision  probability  posed  by  the  disposed  satellite  to  a  secondary  satellite  on  a  circular  orbit  was 
computed  as  a  function  of  the  altitude  of  the  secondary  satellite.  The  collision  probability  is  formulated  from 
the  perspective  of  the  secondary  satellite  as  it  flies  through  the  spatial  density  field  of  the  disposed  satellite. 
The  formulation  for  computing  the  collision  probability  per  unit  time  is  then  derived  by  simplifying  Eq.  1  for 
the  case  when  there  is  no  time  variation  of  the  secondary  satellite  orbit. 

^r-(h)  =  AIcPd(h)v(h)  (5) 

at 

where  h  is  the  altitude  of  the  secondary  satellite,  Pd(h)  is  the  spatial  density  (averaged  over  latitude)  of  the 
disposed  satellite  at  altitude  h,  v(h)  is  the  average  relative  velocity  between  the  disposed  and  secondary 
satellites,  and  Acc  is  the  average  collision  cross-sectional  area  between  the  disposed  vehicle  and  secondary 
satellites.  As  for  the  constellation,  the  spatial  density  Pd(h)  was  computed  using  the  formulation  of 
Dennis/Kessler.  The  value  of  Acc  is  set  to  unity,  so  the  collision  probability  is  normalized  by  collision  cross- 
sectional  area.  The  resulting  formulation  of  collision  probability  assumes  that  the  Q  and  o  for  the  orbits  are 
randomly  and  uniformly  distributed  over  360  deg. 

Figure  5  shows  the  collision  probability  as  a  function  of  secondary  satellite  altitude.  The  eccentricity  of  the 
disposed  satellite  is  varied  parametrically.  The  plot  shows  that,  as  eccentricity  increases,  collision  probability 
is  significantly  reduced  for  most  altitudes  between  the  apogee  and  perigee  of  the  disposed  satellite.  The 
collision  probability  for  an  eccentricity  of  0.36  is  two  orders  of  magnitude  below  the  collision  probability  for 
an  eccentricity  of  0.005  for  most  altitudes.  The  reduction  in  collision  probability  at  the  altitude  extremes  is 
less  substantial,  especially  at  perigee.  This  indicates  that,  for  a  disposal  orbit  undergoing  eccentricity  growth, 
it  is  desirable  to  minimize  the  time  spent  by  the  disposal  orbit  perigee  and  apogee  at  altitudes  with  high 
population  density. 
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Figure  5.  Collision  risk  posed  by  a  disposed  Galileo  satellite  to  a  target  satellite  on  a  circular  orbit  as  a  function  of  target 
satellite  altitude. 

It  is  also  of  interest  to  determine  the  collision  risk  between  two  disposed  satellites  with  the  same  semimajor 
axis,  eccentricity  and  inclination,  but  with  different  and  randomly  oriented  Q  and  00.  For  this  case,  the 
collision  probability  was  formulated  from  the  perspective  of  one  disposed  satellite  as  it  flies  through  the 
spatial  density  field  of  the  other  disposed  satellite.  The  formulation  of  Eqs.  1-3  was  used.  (Eq.  4  was  not 
needed  because  the  in-plane  orbital  elements  of  the  secondary  satellite  are  static  in  this  scenario.)  As  is  done 
throughout  the  study,  the  value  of  Acc  is  set  to  unity. 

Figure  6  shows  the  resulting  collision  probability  as  a  function  of  eccentricity.  The  plot  shows  that,  as 
eccentricity  increases,  collision  probability  between  the  disposed  satellites  is  significantly  reduced.  The 
collision  probability  for  an  eccentricity  of  0.36  is  one  and  a  half  orders  of  magnitude  below  the  collision 
probability  for  an  eccentricity  of  0.005  for  most  altitudes. 


Figure  6.  Collision  risk  between  two  disposed  satellites  as  a  function  of  eccentricity.  Both  satellites  have  the  same 
semimajor  axis,  eccentricity  and  inclination  but  random  Q  and  00. 
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7. 


Disposal  Orbit  Evolution  For  E0  =  0.005 


The  long-term  evolution  was  determined  for  a  disposal  orbit  strategy  with  e0  =  0.005  and  perigee  500  km 
above  the  Galileo  reference  orbit.  The  corresponding  semimajor  axis  is  30647  km.  Initial  inclination  was 
assumed  to  be  56  deg.  The  first  disposal  is  assumed  to  occur  on  1  January  2009.  As  in  the  case  for  the 
constellation  propagation,  Q0  and  cd0  were  varied  over  the  range  [0,360)  deg  in  10  deg  increments.  Orbital 
mean  element  vs.  time  profiles  for  each  (Q0  ,  (Do)  pair  were  generated  by  propagating  over  a  time  period  of 
200  years  using  MEANPROP. 

Figures  7-8  show  the  resulting  evolution  of  apogee  and  perigee  altitude  for  a  sample  set  of  constellation 
orbits  with  Q0  =  0  deg.  Each  pair  of  apogee  and  perigee  curves  corresponds  to  a  specific  value  of  oo0.  From 
this  plot  it  can  be  seen  that  values  of  oa0  are  available  that  result  in  either  large,  moderate,  or  small 
eccentricity  growth.  For  some  of  the  large  eccentricity  growth  cases,  the  perigee  drops  low  enough  so  that  the 
disposed  vehicle  reenters  the  atmosphere  toward  the  end  of  the  200-year  time  period.  In  addition,  for  some  of 
the  large  eccentricity  growth  cases,  the  apogee  altitude  reaches  and  exceeds  the  altitude  of  geosynchronous 
orbit.  However,  a  collision  risk  is  only  posed  to  geosynchronous  satellites  for  relatively  brief  periods  of  time 
when  the  argument  of  perigee  lies  within  one  of  four  very  narrow  bands  that  depend  on  the  current 
eccentricity  of  the  disposal  orbit. 


Figure  7.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.005  for  Q0  =  0  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  8.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  constellation  reference  orbit)  of  Galileo  disposal 
orbits  with  e0  =  0.005  for  Q0  =  0  deg.  Curves  are  shown  for  values  of  oo0  that  span  the  range  [0,360)  deg  in  10  deg 
increments. 


Figures  9-10  show  the  resulting  evolution  of  apogee  and  perigee  altitude  for  a  sample  set  of  constellation 
orbits  with  Q0  =  120  deg.  From  this  plot  it  can  be  seen  that  values  of  030  are  available  that  result  in  either 
moderate  or  small  eccentricity  growth,  but  not  as  large  as  for  the  case  with  Q0  =  0  deg.  There  is  no  case 
where  perigee  drops  low  enough  so  that  the  disposed  vehicle  reenters  within  200  years.  There  is  also  no  case 
where  apogee  reaches  GEO  within  200  years. 
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Figure  9.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.005  for  Q0  =  120  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  10.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  constellation  reference  orbit)  of  Galileo  disposal 
orbits  with  e0  =  0.005  for  Q0  =  120  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg 
increments. 


Figures  11-12  show  the  resulting  evolution  of  apogee  and  perigee  altitude  for  a  sample  set  of  constellation  orbits  with 
Q0  =  170  deg.  From  this  plot  it  can  be  seen  that  values  of  co0  are  available  that  result  in  very  small  eccentricity  growth.  It 
is  not  possible  to  achieve  moderate  or  large  eccentricity  growth. 
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Figure  11.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.005  for  Q0  =  170  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  12.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  constellation  reference  orbit)  of  Galileo  disposal 
orbits  with  e0  =  0.005  for  Q0  =  170  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg 
increments. 
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8. 


The  Minimum  Eccentricity  Growth  Strategy 


Because  of  the  strong  sensitivity  of  eccentricity  growth  to  co0,  a  wide  variety  of  disposal  strategies  are 
potentially  available.  The  strategy  that  has  been  previously  proposed  is  to  minimize  e0  and  select  co0  that 
minimizes  eccentricity  growth.  This  strategy  can  be  called  the  minimum  eccentricity  growth  strategy.  It  will 
minimize  the  number  of  disposed  vehicles  that  penetrate  the  constellation.  By  selecting  the  disposal  orbit 
perigee  to  be  500  km  above  the  constellation  reference  orbit,  the  AV  required  to  clear  the  constellation  is 
minimized. 

There  are  several  disadvantages  with  implementing  this  strategy.  One  disadvantage  is  that  for  small  values  of 
initial  eccentricity,  co0  becomes  very  sensitive  to  small  maneuver  errors,  and  it  is  difficult  to  accurately 
achieve  a  desired  value  of  a)0. 

Another  disadvantage  is  that  disposed  vehicles  will  accumulate  over  time  within  a  tightly  confined 
graveyard.  The  disposed  vehicles  will  experience  in-track  drift  relative  to  each  other  due  to  initial  orbit 
insertion  errors  and  various  orbital  perturbations,  such  as  solar  radiation  pressure.  As  a  result  the  disposed 
satellites  will  pose  a  collision  risk  amongst  themselves. 

In  order  to  model  the  graveyard  population  growth,  a  spatial  density  model  of  the  graveyard  was  constructed 
as  follows.  At  each  value  of  Q0  over  the  range  [0,  360)  deg,  the  value  of  co0  was  selected  that  minimizes  the 
eccentricity  at  the  time  point  100  years  after  disposal.  The  spatial  density  fields  for  each  value  of  Q0  and  each 
time  point  over  the  200-year  period  were  then  computed  and  averaged  together.  The  resulting  averaged 
spatial  density  field  represents  one  disposed  satellite.  To  obtain  the  spatial  density  field  that  accounts  for  all 
satellites  disposed  up  to  a  time  t,  the  single  satellite  spatial  density  field  was  then  multiplied  by  the  number  of 
accumulated  disposed  vehicles  at  time  point  t.  This  model  represents  a  statistically  stationary  model  of  the 
graveyard  in  terms  of  its  orbital  dynamics,  and  is  sufficiently  accurate  for  the  minimum  eccentricity  growth 
disposal  strategy.  It  is  not  accurate  for  strategies  that  result  in  significant  eccentricity  growth,  and  is  not  used 
for  such  strategies  in  this  study. 

Figure  13  shows  the  spatial  density  field  of  the  graveyard  after  100  years,  assuming  a  satellite  disposal  rate  of 
2  satellites  per  year.  It  is  shown  along  with  the  spatial  density  field  of  the  constellation.  It  is  seen  that  the 
graveyard  spatial  density  is  much  higher  than  the  constellation  spatial  density. 
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Figure  13.  Spatial  population  densities  of  the  graveyard  and  constellation  after  100  years  for  the  minimum  eccentricity 
growth  strategy. 

The  cumulative  collision  risk  between  the  disposed  vehicles  in  the  graveyard  was  computed  as  a  function  of 
time.  This  was  done  as  follows.  The  collision  probability  time  profile  was  computed  for  a  disposed  satellite 
that  flies  through  the  spatial  density  field  of  a  previously  disposed  satellite.  The  averaged  single-satellite 
spatial  density  field  that  was  determined  for  the  graveyard  was  used  to  represent  the  previously  disposed 
satellite.  A  collision  probability  profile  was  computed  in  this  manner  for  each  mean  element  vs.  time  profile 
in  the  ensemble  that  was  used  to  generate  the  graveyard  spatial  density  field. 

The  effect  of  a  sequence  of  progressive  satellite  disposals  was  then  modeled.  For  each  subsequently  disposed 
satellite,  a  collision  probability  vs.  time  profile  was  randomly  selected  from  the  ensemble  of  profiles  (thereby 
modeling  the  random  selection  of  Q0)  and  time-shifted  to  start  at  the  time  of  the  current  disposal.  Because 
that  delayed  collision  probability  time  profile  only  accounts  for  the  risk  posed  by  one  previously  disposed 
vehicle,  it  is  first  multiplied  by  the  number  of  previously  disposed  satellites  and  then  added  to  the  previously 
accumulated  collision  probability  time  profiles.  This  process  is  repeated  until  all  vehicles  disposed  over  the 
200-year  period  have  been  included.  The  disposal  sequence  is  then  repeated  1000  times  in  Monte  Carlo 
fashion. 

This  procedure  yields  an  estimate  of  collision  risk  vs.  time  that  is  sufficiently  accurate  for  this  study  because 
the  eccentricity  of  the  disposed  satellites  does  not  grow  significantly.  It  accounts  for  the  effect  of  the 
restricted  dynamical  variation  of  eccentricity  over  time  on  the  graveyard  spatial  density. 

Figure  14  shows  the  cumulative  collision  risk  over  time  between  disposed  satellites  in  the  graveyard.  For 
comparison,  the  collision  probability  between  the  disposed  satellites  and  the  operational  constellation  is  also 
shown.  Both  curves  are  the  median  profiles  resulting  from  the  Monte  Carlo  repetition  of  the  disposal 
sequence.  It  is  seen  that,  for  most  of  the  200-year  time  period,  the  intra-graveyard  collision  risk  is  three  to 
four  orders  of  magnitude  higher  than  the  collision  risk  posed  by  the  disposed  vehicles  to  the  constellation.  As 
an  example  of  the  absolute  collision  risk,  if  the  satellites  have  a  collision  cross-sectional  area  of  100  m2 
(corresponding  to  a  collision  radius  of  5.64  m),  the  probability  that  a  collision  will  occur  sometime  during  the 
200-year  interval  will  be  0.02  (odds  of  one  in  50). 
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Figure  14.  Cumulative  collision  risk  over  time  between  disposed  satellites  in  the  graveyard.  For  comparison,  the 
collision  probability  between  the  disposed  satellites  and  the  operational  constellation  is  also  shown. 


If  a  collision  occurs  between  disposed  vehicles  in  a  graveyard  that  is  located  near  the  constellation,  it  is  very 
likely  that  debris  resulting  from  the  collision  will  follow  orbits  that  penetrate  the  constellation.  This  is  due  to 
the  high  relative  velocities  (typically  on  the  order  of  4.6  km/s)  between  disposed  vehicles  in  combination 
with  the  relatively  weak  gravity  field  in  MEO.  Most  of  this  debris  will  be  smaller  than  the  detection  size 
threshold  of  tracking  networks,  and  therefore  it  will  not  be  possible  for  the  constellation  satellites  to  mitigate 
the  risk  they  pose  via  collision  prediction  and  avoidance.  In  contrast,  the  direct  collision  risk  posed  by  intact 
disposed  vehicles  to  the  constellation  can  be  mitigated  via  collision  prediction  and  avoidance.  Therefore, 
collision  risk  within  the  graveyard  poses  a  significant  indirect  risk  to  the  operational  constellation. 

One  possible  option  to  reduce  the  intra-graveyard  collision  risk  is  to  require  that  subsequently  disposed 
vehicles  move  to  higher  disposal  orbits.  The  disadvantage  of  this  strategy  is  that  it  will  place  a  burden  on  the 
replacement  satellites  that  will  grow  indefinitely  over  time. 
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9.  High  Eccentricity  Growth  Strategies 

The  option  considered  here  is  to  exploit  natural  dilution  of  the  intra-graveyard  collision  risk  by  increasing 
eccentricity  growth.  The  strategy  that  is  investigated  here  is  to  select  (e0,  coo)  to  maximize  eccentricity  growth 
at  100  years  after  disposal  or,  if  atmospheric  reentry  can  be  achieved,  to  minimize  time  to  reentry.  The 
disposal  orbit  perigee  is  anchored  at  500  km  above  the  constellation  reference  orbit.  An  advantage  of  this 
strategy  is  that,  as  e0  increases,  it  becomes  easier  to  accurately  target  co0-  A  disadvantage  is  that  increasing  e0 
must  be  accomplished  by  raising  apogee.  This  will  cost  extra  AV,  so  more  effort  on  the  part  of  the  satellite 
will  be  required  in  the  form  of  propellant  as  e0  is  increased.  An  important  issue  to  consider  is  that  the  high 
eccentricity  growth  strategy  will  pose  a  higher  direct  collision  risk  to  the  constellation  than  the  minimum 
eccentricity  growth  strategy.  Therefore,  the  tradeoff  between  graveyard  collision  risk  and  constellation 
collision  risk  must  be  considered. 

The  high  eccentricity  growth  strategy  will  also  pose  a  higher  direct  collision  risk  to  other  objects  external  to 
Galileo,  such  as  the  GPS  and  Glonass  systems  and,  in  extreme  cases,  low  Earth  orbit  and  geosynchronous 
objects.  Determination  of  the  collision  risk  posed  to  these  objects  was  beyond  the  scope  of  this  study. 
However,  it  is  expected  to  be  low  due  to  the  dilutional  effect  of  eccentricity  growth  as  illustrated  in  Figs.  5-6. 
It  is  planned  to  quantify  the  comprehensive  collision  risk  between  all  disposed,  operational,  and  background 
objects  in  a  future  study. 

Four  levels  of  eccentricity  were  considered.  Table  1  shows  the  corresponding  disposal  orbit  semimajor  axis 
and  the  required  increase  in  semimajor  axis  above  the  constellation  reference  orbit  for  each  level  of 
eccentricity.  For  reference,  e0  =  0.021  will  have  an  increase  in  semimajor  axis  of  1156  km.  In  comparison, 
the  increase  in  semimajor  axis  for  historical  GPS  disposals  has  ranged  from  750  to  1350  km,  with  a  range 
midpoint  of  1050  km  (Jenkin  and  Gick,  2003).  Therefore,  e0  =  0.021  will  have  an  effort  that  is  roughly 
comparable  to  that  of  the  mid-range  of  historical  practice  of  GPS  disposal.  Eccentricities  e0  =  0.036593  and 
0.051649  will  require  roughly  1.5  to  2  times  the  effort,  respectively,  of  the  mid-range  of  historical  GPS 
practice.  This  range  of  disposal  effort  was  selected  for  this  study  because  it  is  believed  to  be  feasible  with 
modem  satellite  technology. 


Table  1.  High  eccentricity  growth  strategy  cases. 


eO 

Perigee  altitude 
rel.  to  Galileo 
(km) 

Semimajor 
axis  (km) 

Increase  in 
semimajor  axis 
above  Galileo  (km) 

0.005 

500 

30647 

653 

0.021051 

500 

31150 

1156 

0.036593 

500 

31652 

1658 

0.051649 

500 

32155 

2161 

The  resulting  apogee  and  perigee  altitude  profiles  for  e0  =  0.005  have  already  been  presented  in  Figs.  6-11. 
As  noted  before,  large  eccentricity  growth,  as  well  as  re-entry  within  200  years,  can  be  achieved  for  some  but 
not  all  values  of  Q0- 
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Figures  15-17  show  the  apogee  and  perigee  altitude  profiles  for  e0  =  0.021  and  Q0  =  0,  120,  and  170  deg 
respectively.  It  is  seen  that  eccentricity  growth  has  been  increased,  even  for  Q0  =  170  deg.  In  addition,  for  Q0 
=  0,  re-entry  can  be  achieved  within  140  years  after  disposal,  and  for  Q0  =  120,  re-entry  can  be  achieved 
within  200  years  after  disposal. 


Figure  15.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.021  for  Q0  =  0  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  16.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.021  for  Q0  =  120  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  17.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.021  for  Q0  =  170  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 


Figures  18-20  show  the  apogee  and  perigee  altitude  profiles  for  e0  =  0.036593  and  Q0  =  290,  120,  and  170 
deg  respectively.  It  is  seen  that  eccentricity  growth  has  been  increased  further,  especially  for  Q0  =  170  deg. 
For  Q0  =  120  deg,  re-entry  can  be  achieved  within  180  to  190  years  after  disposal.  For  Q0  =  290  deg,  re-entry 
can  be  achieved  as  early  as  120  years  after  disposal. 
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Figure  18.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.036593  for  Q0  =  290  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  19.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.036593  for  Q0  =  120  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  20.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.036593  for  Q0  =  120  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 


Figures  21-23  show  the  apogee  and  perigee  altitude  profiles  for  e0  =  0.051649  and  Q0  =  0,  120,  and  170  deg 
respectively.  It  is  seen  that,  for  Q0  =  0  deg,  the  perigee  altitude  profiles  undergo  a  reversal,  and  re-entry  does 
not  occur  within  200  years.  However,  for  Q0  =  120  deg,  re-entry  can  be  achieved  130  years  after  disposal. 
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Figure  21.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.051649  for  Q0  =  0  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  22.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.051649  for  Q0  =  120  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 
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Figure  23.  Apogee  and  perigee  altitude  profiles  (shown  relative  to  the  Earth  surface)  of  Galileo  disposal  orbits  with  e0  = 
0.051649  for  Q0  =  170  deg.  Curves  are  shown  for  values  of  co0  that  span  the  range  [0,360)  deg  in  10  deg  increments. 


Figure  24  shows  the  time  profiles  of  collision  risk  posed  to  the  constellation  by  accumulating  disposed 
vehicles  (assuming  a  disposal  rate  of  two  per  year)  corresponding  to  the  various  disposal  strategies.  The  plot 
shows  curves  for  the  four  high  eccentricity  growth  cases  considered,  and  also  a  strategy  in  which  e0  =  0.005, 
perigee  is  500  km  above  the  constellation  reference  orbit,  and  oo0  is  not  specifically  targeted  but  rather 
allowed  to  vary  randomly  among  sequentially  disposed  satellites.  From  the  plot,  it  is  seen  that  the  case  with 
the  highest  collision  risk  is  the  high  eccentricity  growth  case  with  e0  =  0.005.  The  case  with  the  second 
highest  collision  risk  after  130  years  is  the  random  cd0  case  with  e0  =  0.005.  The  long-term  collision  risk  is 
reduced  further  for  the  high  eccentricity  growth  strategy  cases  with  higher  e0.  The  collision  risk  is  continually 
reduced  as  e0  is  increased,  however  the  gains  in  risk  reduction  begin  to  diminish  after  e0  exceeds  0.036593. 
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Figure  24.  Collision  risk  posed  to  the  Galileo  constellation  by  various  disposal  strategies. 
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Figure  25  shows  the  collision  risk  posed  to  the  constellation  for  the  disposal  strategies  considered,  including 
the  minimum  eccentricity  growth  case,  along  with  the  intra-graveyard  collision  risk  for  the  minimum 
eccentricity  growth  case.  It  is  necessary  to  use  a  logarithmic  plot  in  order  to  show  the  broad  range  of 
collision  risk  that  is  encompassed.  It  is  seen  that  the  high  eccentricity  growth  strategies  and  the  random  co0 
strategy  result  in  a  higher  collision  risk  to  the  constellation  than  the  minimum  eccentricity  growth  strategy. 
However,  the  collision  risk  posed  to  the  constellation  by  the  high  eccentricity  growth  strategies  is  much 
lower  than  the  intra-graveyard  collision  risk.  The  high  eccentricity  growth  strategy  with  e0  =  0.036593 
produces  a  constellation  risk  that  is  more  than  two  orders  of  magnitude  below  the  graveyard  risk  for  the 
minimum  eccentricity  growth  strategy. 


Figure  25.  Comparison  of  the  collision  risk  posed  to  the  constellation  by  various  disposal  orbit  strategies  with  the  intra¬ 
graveyard  collision  risk  for  the  minimum  eccentricity  growth  strategy. 
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The  intra-graveyard  collision  risk  for  the  high  eccentricity  growth  strategies  was  not  quantified  in  this  study. 
However,  from  Fig.  6  it  is  expected  that  these  strategies  will  significantly  reduce  the  intra-graveyard  collision 
risk.  Therefore,  the  combined  constellation  and  graveyard  collision  risk  associated  with  Galileo  can  be 
significantly  reduced  by  using  the  high  eccentricity  growth  strategy  as  long  as  e0  is  large  enough. 

Table  2  shows  the  percentage  of  disposed  satellites  that  will  re-enter  within  200  years  after  disposal  for 
various  strategies.  The  results  account  for  the  entire  [0,360)  deg  range  of  Q0-  It  is  seen  that  the  re-entry 
percentage  is  significantly  increased  and  is  substantial  for  the  high  eccentricity  growth  cases  with  e0  =  0.021 1 
and  larger. 


Table  2.  Percentage  of  satellites  that  re-enter  within  200  years  after  disposal. 


(D0-selection  strategy 

e0 

%  of  satellites  that  reenter 
before  200  years  after 
disposal 

Random  co0 

0.005 

7.1% 

High  eccentricity 

0.005 

19% 

IV 

0.021 1 

67% 

IV 

0.036593 

75% 

M 

0.051649 

92% 
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10. 


Conclusions 


The  study  presented  here  investigated  the  possibility  of  diluting  MEO  disposal  orbit  collision  risk  by 
selection  of  initial  disposal  orbit  eccentricity  and  argument  of  perigee  in  order  to  control  long-term 
eccentricity  growth.  The  Galileo  constellation  was  selected  as  an  example. 

Various  disposal  strategies  were  considered.  The  minimum  eccentricity  growth  strategy  minimizes  the 
collision  risk  posed  directly  to  the  constellation  by  accumulating  disposed  vehicles.  However,  it  also  results 
in  a  risk  of  collision  between  disposed  vehicles  in  the  graveyard  that  is  much  higher  than  the  risk  posed 
directly  to  the  constellation.  The  minimum  eccentricity  growth  strategy  also  has  the  disadvantage  that  it  is 
difficult  to  accurately  target  argument  of  perigee,  thereby  making  it  difficult  to  attain  the  desired  eccentricity 
growth  profile. 

High  eccentricity  growth  strategies  increase  the  collision  risk  posed  directly  to  the  constellation,  but  they  will 
significantly  reduce  the  intra-graveyard  collision  risk.  Therefore,  the  high  eccentricity  growth  strategies  will 
reduce  the  combined  constellation  and  intra-graveyard  collision  risk.  The  increase  in  collision  risk  posed  to 
the  constellation  can  be  reversed  by  increasing  initial  disposal  orbit  eccentricity.  The  high  eccentricity 
growth  strategies  also  have  the  advantage  that  it  is  easier  to  accurately  target  argument  of  perigee,  thereby 
making  it  easier  to  attain  the  desired  eccentricity  growth  profile. 

It  is  preferable  to  avoid  a  high  graveyard  collision  risk  because  potential  graveyard  collisions  can  produce 
untrackable  debris  that  will  penetrate  the  constellation  and  cannot  be  evaded  by  collision  avoidance. 

Determination  of  the  collision  risk  posed  to  objects  external  to  the  Galileo  system  was  beyond  the  scope  of 
this  study.  However,  it  is  expected  to  be  low  due  to  the  dilutional  effect  of  eccentricity  growth.  It  is  planned 
to  quantify  the  comprehensive  collision  risk  between  all  disposed,  operational,  and  background  objects  in  a 
future  study. 

High  eccentricity  growth  strategies  also  offer  the  option  of  significantly  increasing  the  percentage  of 
disposed  vehicles  that  will  re-enter  the  atmosphere  within  200  years  after  disposal  rather  than  remain  on  orbit 
for  thousands  of  years.  This  option  may  be  desirable  if  the  vehicles  are  designed  to  pose  an  acceptably  low 
ground  casualty  expectation. 

Therefore,  both  in  terms  of  collision  risk  reduction  and  removal  of  disposed  vehicles  from  orbit,  the  strategy 
of  achieving  high  eccentricity  growth  by  selection  of  initial  disposal  orbit  eccentricity  and  argument  of 
perigee  offers  an  effective  and  potentially  inexpensive  option  for  MEO  debris  mitigation. 
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